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In the presence of CP violation, the effective Hamiltonian matrix describing a neutral meson 
anti-meson system does not commute with its hermitian conjugate. As a result, this matrix cannot 
be diagonalized by a unitary transformation and one needs to introduce a reciprocal basis. Although 
known, this fact is seldom discussed and almost never used. Here, we use this concept to highlight 
a parametrization of the Hamiltonian matrix in terms of physical observables, and we show that 
using it reduces a number of long and tedious derivations into simple matrix multiplications. These 
results have a straightforward application for propagation in matter. We also comment on the 
(mathematical) relation with neutrino oscillations. 
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I. INTRODUCTION 

We are interested in the efFective 2x2 Hamiltonian matrix describing the mixing in the P° 
P stands for K, D, B^, or Bg. We denote this 2x2 matrix hy H ^ M ~ i/2T where 



M 



H + W) /2 and - ir/2 ^ [H - W ] /2, 



systems, where 



(1) 



describe the hermitian and anti-hcrmitian parts of H , respectively. Both M and F are hermitian. Matrices satisfying 
are called "normal" matrices. It is easy to show that 



t _ 



= ^ [M,r] = 0. 



(2) 



Moreover, a matrix is normal if and only if it can be diagonalized by a unitary transformation. It is often stated that 
non-unitary transformations arise whenever H is not hermitian. This is not the case. What is relevant is whether H 
is normal or not. Indeed, if F 7^ then H is not hermitian; however, H can still be diagonalized by a unitary matrix 
as long as [M^, F] = 0. 

In section ^ we introduce the concept of 'reciprocal basis' and we show that the presence of T violation in the 
po _ po system forces us to use such a basis. The physical observables are defined in section III and they are used 
in section to parametrize H exclusively in terms of measurable quantities. The time evolution of the P° — 
system is di scuss ed in section Section VI explains why the P° — P*^ should be considered as intermediate states, 
and section |VII| sho ws an error which arises when one does not use the reciprocal basis. Matter effects are then 
considered in section VII] . This differs from all previous analyses of matter effects in that no use is made of the Good 
equations; here the time evolution is obtained in a trivial way. In section [X we compare the mixing in the P^ — P^ 



system with the mixing in the neutrino sector. To this end, we start by showing how the equation describing the 
time evolution and its solution would look if we had chosen a different reference frame. We present our conclusions in 
section For completeness, appendix ^ contains some elementary notions of collision theory, which are needed to 
describe the evolution in the presence of interactions with matter. Appendix ^ contains two other parametrizations 
of the physical observables commonly found in the literature, the first of which is most convenient for the comparison 
with the neutrino sector. 



II. THE RECIPROCAL BASIS 



*Permanent address: Instituto Superior de Engenharia de Lisboa, Rua Conselheiro Emidio Navarro, 1900 Lisboa, Portugal. 



1 



A. Definition 



Why do we change basis at all? One reason is that the time evolution of the state \'4'{t)) describing the — 
mixed state, which is given by 

ij^\m) = H\m). (3) 

becomes trivial in the basis in which H is diagonal. Eq. and H have been written in the P° — P° rest frame and 
t is the proper time. 

We denote the (complex) eigenvalues of H hy fia — Tna — i/'2^a and /ih = mf, — i/2Tiy, corresponding to the 
eigenvectors 

\Pa)\_(pa qa \ ( ^ \ _ ( ^ \ 

As a result, the matrix H is diagonalized through 

X^HX =( ^ V (5) 



fib 



where 



X-' = ). (6) 

PaQb + PbQa \ la -Pa J 

As stated above, H is normal if and only if X is unitary. This is what one learns in algebra. 

So, why do (most) people worry about performing non-unitary transformations? The reason is that one would like 
the mass basis {iJ'a), \Pb)} to retain a number of the nice (orthogonality) features of the {l-P*^), \P^)} flavor basis. 
Among these: the orthogonality conditions 

(pO|P") = (P0|P°) =0, 

(P°|pO) = (P0|p0) = 1; (7) 
the fact that |P*')(P''| and |P'^)(P°| are projection operators; the completeness relation 

|P0)(P0| + |P^)(P0| = 1; (8) 
and the decomposition of the effective Hamiltonian as 

n = |P°)i/ll(P°| + |P°)i/l2(n + |P^)i?2l(P"| + \P^)H22{P^\ 

= (|P0), |po))jj(^(^j^. (9) 

All these relations involve the basis of flavor eigenkets { | P°) , | P°) } and the basis of the corresponding bras { (P'' | , (P"^ | } . 
The problem is that, when H is not normal, we cannot find similar relations involving the basis of mass eigenkets 
{I Pa), |Pb)} and the basis of the corresponding bras, {(Pal, (A|}- In particular, it is easy to see from the diagonalization 
in Eq. (||) that the analogue of Eq. (||) is 



n=\Pa)fla{Pa\ + \Pb)fib{Pb\ 
fJ-a 

Aifc ; V {Pb\ 



Pa)^ \p») ) ( ,M ( ] ■ (10) 



This does not involve the bras (Pal and (Pb|, 



2 



but rather the so called 'reciprocal basis' 



The reciprocal basis may also be defined by the orthogonality conditions 

(Pain) = {Pb\Pa) = 0, 

{Pa\Pa) - {Pb\Pb) = 1. (13) 

Moreover, \Pa){Pa\ and |Ph)(P{,| are projection operators, and the partition of unity becomes 

\Pa){Pa\ + \Pb){Pb\^l. (14) 

If H is not normal, then X is not unitary, and {(-Pal, {Pb\} in Eq. ( pi] ) do not coincide with |(Pa|, (Alj in Eq. (p^). 

Another way to state this fact is to note that H is normal {X is unitary) if and only if its right-eigenvectors coincide 

with its left-eigenvectors. 

That these features have an impact on the K'^ — K° system, was pointed out long ago by Sachs and by Enz 
and Lewis [^. More recently, they have been stressed by Alvarez-Gaume et al. Q|, and by Branco, Lavoura and Silva 
in their book "CP violation" js). Still, we have found that they are not common knowledge. This is unfortunate 
since there are a number of results that usually require considerable algebra which become trivial once the matrix 
formulation discussed here is implemented. Moreover, one can express the matrix elements of iJ, written in the 
po _ po basis, in terms of observable quantities. This is what we show here. 



B. The relation to CP violation 

We will now show that the reciprocal basis is required by the observation of T and CP violation in the mixing in 
the neutral meson systems. The discrete symmetries have the following effects on the matrix elements of H: 

CVT conservation =^ Hn = H22, 
T conservation^ 1-^121 = 1-^21 1, 
CV conservation => Hu = H22 and \Hi2\ = \H2i\- (15) 

The 1964 discovery that \Hi2\ 7^ 1-^21 1 in the kaon system Q means that there is T and CP violation in K'^ — K'-' 
mixing. Moreover, since the (1,1) entry in the matrix [H,H''] is given by |-ffi2p — 1-^21^, this experimental result 
also implies that the matrix H is not normal and, thus, that we are forced to deal with non-unitary matrices in the 
neutral kaon system. 

For the other neutral meson systems, |i?i2| 7^ 1-^21 1 has not been established experimentally. Nevertheless, the 
Standard Model predicts that, albeit the difference is small, |-ffi2| 7^ 1-^211 does indeed hold. As before, this implies 
CP violation in the mixing and forces the use of the reciprocal basis in all the neutral meson systems. 

III. OBSERVABLES IN THE pO po MIXING 

Let us start by introducing some notation. We define 

fi = m- iT/2 = {fia + f^b)/'^, 
= Am - iAr/2 = /ia - /Zfc. (16) 

Sometimes it is convenient to trade the eigenvalue difference for x — iy = A^/T. We may write the mixing matrix X 
in terms of new parameters 

9a _ gb 
n _ Pa Pb 



Pa Pb 

and 



3 



q 
p 



PaPb 



(18) 



Notice that we have not defined the quantities q and p separately; we only define the ratio q/p. With this notation 
the mixing matrix may be re-written as 



X = 




1 



Pa 

1-0 pV i+e I \ ^ 



(19) 




(20) 



We point out that these transformation matrices involve the normalization constants Pa and pb- Finally, it will also 
prove convenient to define 



(21) 







2 


1 - 


9 






P 








2 


1 + 


9 






P 





meaning that \q/p\ = \J\^- 

The fact that the trace and determinant are invariant under the general similarity transformation in Eq. (^ implies 
that 



Moreover, from 



A^i = ^AHi2H2i + {H22 - HiiY 



(22) 



we find that 



leading to 



/ Hii H12 




\ H21 H22 




Hii H12 \ 


f Pb 


H21 H22 J 


\ -qb 


qa P-a 





fJ'a 



t^b 



Pa 

qa 

Pb 

-qb 



21 



Pa 

qb_ 
Pb 



H12 

Hii — fib 



H 



12 



/ia — H22 ' 

H2I 
H22 — Mb ' 



22 



fJ.a ' 
\Hl2\ 



fJ^b 



\H2l\ 



121 



H 



211 



(23) 



(24) 



(25) 



and q/p = \/H2i/Hr2. We see that Re 6* and Im6' arc CP- and CPT-violating, while 5 is CP- and T- violating. 

Although H contains eight real numbers, only seven are physically meaningful. Indeed, one is free to change the 
phase of the kets |pO), |P0), |Pq), and |Ph), as 



\p')- 


^ e*'^|P°), 




^ e*~|PO), 


\Pa)^ 




in)- 


-^e^^^lP,) 



(26) 
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Under these transformations 



q/p^e^^-<-'^\/p, (27) 

while Hii, H22, /i, A/i, 6*, and 5 do not change. Therefore, the relative phase between H12 and H21 is physically 
meaningless and H contains only seven observables. Similarly, the phase of q/p is also unphysical. As a result, we 
have four observables in the eigenvalues, p and A^, and three in the mixing matrix, 9 and S (or, alternatively, \q/p\)- 



IV. PARAMETRIZING H WITH MEASURABLE QUANTITIES 

Eqs. ( p2|) and ( [25| ) give the measurable mixing and eigenvalue parameters in terms of the Hij matrix elements 
which one can calculate in a given model. Given the current and upcoming experimental probes of the various neutral 
meson systems, it seems much more appropriate to do precisely the opposite; that is, to give the Hij matrix elements 
in terms of the experimentally accessible quantities. Such expressions would give My and Tij in a completely model 
independent way, with absolutely no assumptions. One could then calculate these quantities in any given model; if 
they fit in the allowed ranges the model would be viable. 

Surprisingly, this is not is done in most expositions of the P'^ — P'^ mixing. The reason is simple. Eqs. ( [2^ ) 
and (|2^) are non-linear in the Hij matrix elements. Thus, inverting them by brute force would entail a tedious 
calculation. With the matrix manipulation discussed here this inversion is straightforward. Indeed, Eq. (|^) can be 
trivially transformed into |^] 



H = X 



/ia 





(28) 



where we have used Eqs. (|T^) and (pO|). This equation expresses in a very compact form the relation between the 
quantities which are experimentally accessible and those which are easily calculated in a given theory. Expanding it, 
we find 

Am Ar 

Mil = m - Rc6' Im6' , 

2 4 



Mo 



-Ml 



Am Ar 

m + Re 6* h Im6' , 

2 4 



1 



2(1 + <5) 



Re(Vl" 



Am — iS ■ 



Ar 



Ar 



+ iS Am 



(29) 



and 



Ar 

rii = r - Re0 — + Im6'Am, 



r22 = r + Re( 



Ar 



Im0 Am, 



'^r - ^ 
-J- 12 — 

p 



1 



Ar 



Re(vT^^) ( — +«(^Amj -Im(Vr^^) f Am-i5 ^ 



Ar 



(30) 



We would argue that this is the best way to quote the experimental results. The impact of any assumption made 
about the physical observables, such as CPT or T conservation, is transparent in Eqs. (|29| ) and (^0[). 

A few remarks are in order. Firstly we note that Eqs. (|l^ ) and ( ^o|) involved the overall normalization factors Pa 
andpb, but that these cancel in the multiplication on th e rig ht hand side of Eq. (p8|). Secondly, although M12, ri2 and 
q/p are not rephasing invariant, we can see from Eqs. (27) that q/pMi2, q/pTi2 and Afi2 1^2 indeed physically 
meaningful. Thirdly, the equations involving T are needed also for the unitarity conditions 



5 



^I^gl^ =rii =r(l-yRe6' + xIm( 

g 

^141' =r22 =r(l + yRe0-xIm( 



a 



^« - q ^ {y + iSx) RejVr^) - {x - iSy) ImiVT^) 



g 



where Ag = (g|T|P°), Ag — {g\T\P^), and the sums run over aU the available decay modes g. 

V. TIME EVOLUTION 

The time evolution of the neutral meson system is easily obtained using Eqs. (|lO|) and and the fact that 

\Pa){Pa\ and \Pb){Pb\ are projection operators. We find, 

expi-iHt) = e-'''^'\Pa){Pa\+e-'^'*\Pb){Pb\. 



It is now trivial to write the evolution operator back in the flavor basis. Indeed, using Eqs. (jj) and (|12|), we find 
exp i-^m) = ( |P"), \n)xi^ J^^, ) j^j ) 



where 



1 r cos(^) 

g±{t) EE - (e-^''"* ± e-'""') = e"*™* e'^^^ <^ 1 ' (^4) 

2 [ — isin (— f-^) 

This corresponds to the usual expressions for the time evolution of a state which starts out as P° or P^ 

\P°{t)) = exp {-^m)\P') = [g+it) 9g^{t)] |P°) + ^VT^9-{t) \n, 

P 

\P^{t)) = exp (-iWi)lPO) = P^^/r~¥g_{t)\P°) + [g+{t) + 0g^{t)] |P0), (35) 

respectively. At this point it is important to emphasize the fact that, in deriving this result, no assumptions were 
made about the form of the original matrix H . This observation will become important once we consider the evolution 
in matter. 

VI. NEUTRAL MESONS AS INTERMEDIATE STATES 

Because there is CP violation in P° — P" mixing, there is no selection rule allowing us to choose a final state / to 
which Pa (or P^) can decay while Pf, (P^) cannot. That is, all calculations must involve the full transition chain |^ 

i^X{Pa,Pb}-^Xf, (36) 

with both neutral meson eigenstates as intermediate states, in order to be formally correct. Obviously, one could 



ignore this problem. Still, as we show in section VI], one will be lead into incorrect results if the reciprocal basis is 
not used as the 'out' bra. 
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Recently, Amorim, Santos, and Silva JTo| have highlighted a very important point about the transition chain in 
Eq. (|3^). They showed that this evolution can be fully parametrized by the usual quantities A/ and \j, describing 
the decays {P^,P^} — > /, /, supplemented by two new quantities and describing the production mechanism 
i,i ^ {P",pO}. (Although they applied these results only to the case in which i,i {P",pO} represents a decay, 
their formalism is valid in all generality.) The new quantities and may entail new sources of CP violation, just 
like A/ and Ajr do. They are absent from the decays B J/ifiK — > J/^/'[/]/f studied previously 0| because, in those 
cases, the initial B° meson can only decay to one of the kaon's flavor eigenstates. However, they are crucial for the 
decays — D + [f]o + and, in general, whenever the initial state i can produce (or, in particular, 

decay into) both flavor eigenstates of the intermediate neutral meson system, P'^ and P^. 

Let us consider the decay chain i X{Pa,Pb} Xf. The complete amplitude for this process involves the 
amplitude for the initial decay into XPa or XPh, the time-evolution amplitude for this state, given by Eq. (|3^), and 
finally the amplitude for the decay into Xf. Suppressing the reference to X, we find 

A{i^ Pa,b ^ /) = {f\T\Pa) e-"^"' {Pa\m + (/|T|Pfc) e"'''"* {Pb\m. (37) 

This is an exact expression. However, sometimes it is possible to choose a final state / and to set the experimental 
conditions in such a way as to maximize the importance of i ^ XPa X f relative to i — > XPh ^ X f. In that case 
we may make the approximation 

A{i^Pa.b^f) ^A{i^Pa^f) 

= (/|r|p,)e-'^«*(p,|r|z) 

= (/iriPa)e-^^"* [(pjp")(p"iT|z) + {Paininm] , 

(38) 

where we have used the partition of unity |P'')(P*'| + |P")(P°I = 1 to derive the last line. When one uses the 
approximation in Eq. (^), one talks about 'the decay i — s- XPq',Q and writes 

Ait^ XPa) - (Pa|P°) A{t ^ XP") + (Pa|P^> Aii ^ XPO) 

= I [p'^A{i ^ XP°) + q-^A{i ^ XPO)j ^ (39) 

where, in the last line, we have assumed the CPT-invariant case: 

{Pa\ = l{p-\P"\+q-'{n 

{Pb\^l{p-\P'\-q-'{n). (40) 

Therefore, the ratio of the two component amplitudes in Eq. (|39| ) is given by q^^ /p^^ = p/q, and not by q* /p* — as 
would have been the case if we had used {Ph\ instead of {Ph\- The difference between q~^/p~^ and q*/p* only 
disappears in the limit \q/p\ — 1. We will now show that this has a formal impact in the study of the decay 
Bd -> J/ijKs. 



VII. ON THE NEED FOR THE RECIPROCAL BASIS IN Bd ^ J/t/^Ks 

This decay is so important that it is surprising how many times it is calculated without even mentioning that the 
use of the reciprocal basis is required in order to obtain the exact result. We repeat, in this decay the use of the 
reciprocal basis is not a convenient calculational tool. It is unavoidable when one wishes to obtain the result without 
approximations. 



^ Nevertheless, strictly speaking, it is Eq. ( ply which expresses the correct way to think about decays into neutral-meson 
eigenstates [^,^. As we stressed above, the point is that, since CP is violated, there is no final state / that can be obtained 
only from Pa and not from Pb- There will always be a non-zero amplitude for the decay path i XPb Xf. 
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The first observation we should make is that what one looks for experimentally is the decay chain Bd J/ tpK — > 
J/iP{'!T'!t)k, and that both intermediate Ks and contribute to this decay. The following argument should make 
it clear why the intermediate must contribute. Consider the decay chain Bd J/ipK J/'0(7r7r)x, but where 
we have chosen to look only for kaons which live a proper time r ^ T5 before they decay. Clearly, for these 
kaons, the Kg component will have disappeared before the decay, and all tttt final states must have come from an 
intermediate K^. This explains why, in general, one must use Eq. (^7|). However, in the experiments searching 
for Bd — > J/^Ks one is looking at kaon proper times r < IOT5. Therefore, in these experiments the decay path 
Bd J/i/jKl J / 'ip{'K'K) K is very suppressed with respect to the decay path Bd J/^Ks J/'iP(tt'k)k both due 
to the huge ratio A{Ks tttt) / A{K l tttt) and to the time interval probed. This leads us to Eq. ( p8| ) and, ignoring 
the normalizations Pa and ph, allows us to talk about the decay Bd J/ipKs as in Eq. (^9|). 

Having established under which circumstances we may (to good approximation) talk about the decay Bd —>■ J/^Kg, 
we are now in position to describe the upcoming measurement of CP violation in this decay. These experiments will 
determine the imaginary part of 



A 



_ QBd AjBd ^ J/^pKs) 
PBd A{Bd ^ J/i^Ks) ' 



(41) 



We wish to calculate A{Bd J/tpKs) and A{Bd J/tpKs)- We recall that the decays Bd — > J/^pK'^ and Bd — > 
J/tpK'^ are forbidden to leading order in the SM, and, to simplify the problem, we consider the CPT-conserving case, 
in which 

\Ks)=PK\K")-qK\K^), 
{Ks\=p},{K°\~q},m, 

{Ks\^l[pK\K'\-q^'(K^\\. (42) 

The question is whether one should use {Ks\ or {Ks\ in the final state. That is, we wish to know whether to use 

J/^if°) + {ks\K^) A{Bd J/ifK^) 



A {Bd -> J/i^Ks) = (i^sl A'°) A{B, 
_ 1 

~ 2 
1 



p-K^A{Bd J/i/^X°) - q^^AiBd ^ J/ijK^) 



and 



-pl^AiBd ^ J/^if"), 



A [Bd ^ J/^Ks) = -]^q^^A[Bd ^ J/i>K^) 



(43) 
(44) 



or, alternatively, use 



A {Bd J/i^Ks) = {Ks\K°) A{Bd ^ J/V^if") + {Ks\K^) A{Bd ^ J/^jK^) 
= p*KA{Bd -> J/^if°) - q*KA{Bd -> J/^K^) 
^p*KA{Bd^J/^K°), 



and 



A (Bd ^ J/^jKs) = -qKA{Bd -> J/iPK^). 



In the first case we obtain 



in the second we obtain 



A 



qBd AjBd ^ J/i^K^) PK 
PBd A{Bd^ Jl^^K^)qK 



A 



qBd A{Bd ^ J/^RO) 

PBd A{Bd^ J/i^KO)p*K 



(45) 



(46) 



(47) 



(48) 



From the previous section, we know that the first expression is the correct one. And, in deriving it, we had to know 
what the reciprocal basis was and that it had to be used. Nevertheless, since \qK/PK\ only differs from one at order 
lO"'^ and we are looking for a large effect in Asd^j/^z-Xs: this detail, although needed for an exact formal derivation, 
is numerically insignificant. This explains why it has gone largely unnoticed |13[. 
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VIII. MATTER EFFECTS IN THE P° po EVOLUTION 

We now wish to study how the time evolution of the P° — changes in the presence of matter. It should be clear 
that the matter effects will change the specific form of H but, since we have considered the most general such matrix, 
all the derivations presented above should still apply. It remains to relate the parameters in matter and in vacuum. 

We will denote the matrices, matrix elements and eigenvalues in vacuum by unprimed quantities and their analogue 
in matter by primed quantities. For example, when kaons transverse matter, they are subject to strong interactions 
which conserve strangeness but which treat the K'^ and differently.^ This effect may be parametrized by a new 
effective Hamiltonian 

^nuc=(j^), (49) 

which must be added to the Hamiltonian in vacuum. Notice that this parametrization is completely general. It 
describes any strangeness-conserving interaction whatsoever. It is also important to notice that our original evolution 
equation, Eq. (|^), and vacuum Hamiltonian H have been written in the — P° rest frame. Before we add i?nuc to 
H we must ensure that H^uc is also expressed in the rest frame. This point is discussed in appendix 
The full Hamiltonian in matter becomes 

H' = H + H^^,. (50) 



Now, we have already studied the most general effective Hamiltonian, and Eq. ( pSf ) relates such an Hamiltonian written 
in the flavor basis with the corresponding eigenvalues and mixing parameters. Therefore, relating the observables in 



vacuum and in matter becomes another simple exercise. Eqs. (28), ([49| ) and ( pOD yield 



(51) 



A few features are worth mentioning. Firstly, = -^12 &nd iJji = H21. As a result, q' /p' — q/p- In particular, 
the CP- and T-violating parameter 5, which depends on \q' /p'\ = \q/p\: is the same in vacuum and in the presence of 
matter. Therefore, the parameters in vacuum and in matter are related through. 



/ , X + X 



A/i' = V(Am)^ + 26* A/i Ax + (Ax)2 = A/i y/T+irT+ir^, 

^, _ A^6l + Ax e + 2r 

" ^(A/i)2 + 261 A^ Ax + (Ax)2 ^ Vl + 4r0 + 4r2' 



(52) 



where Ax = X~Xy ^iid we have introduced the 'regeneration parameter' r = Ax/(2A/i). It will also prove convenient 
to find 



^(A^)2 + 2e'A^Ax + (Ax)2 V l + 4r6'-f 4r2- 



Secondly, it is clear from Eq. (51), and also from Eqs. ( p2| ) and (|53|), that the flavor-diagonal matter effects considered 
here act just like violations of CPT. Thirdly, we expect the matter effects to be much larger than any (necessarily 
small) CPT-violation that there might be already present in vacuum. Therefore, we may take 6* = to get 

/ , X + X 
= A^H 7^, 



A^i' = v/(AAi)2 + (x-x)' = A/i Vl + 4r2, 
0' ^ X- X ^ 2r 

v/(Aa.)2 + (x-x)2 VTT4^' 



(54) 



^The total cross-section for K'' interacting with a nucleus is larger than that for JC" on the same nucleus. For example, 
K°p Att"*" takes place but there is no corresponding reaction for K'\ 
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and y/1 — 0'"^ = l/Vl + 4r^. We stress that Eq. ( ^l| ) is completely general, as will be the time evolution based on it. 

The time evolution in matter is now trivial to find. It is given in Eqs. [or, alternatively, in Eqs. (^] and (HJ), 
with the unprimed quantities substituted by the primed quantities. This solution had been found for the kaon system 
by Good jlj], building on earlier work by Case ^5|. Recent re-derivations may be foimd in references and 
In all these articles, the authors write a new evolution equation obtained by combining the diagonalized form of H 
with the new term i?nuc written in the {Kl,Ks} basis. Thus, they would seem to be solving a new complicated 
set of equations: the so-called 'Good equations'. In the method presented here, we have made no reference to 'new' 
differential equations. We had already solved the most general evolution equation once and for all, Eqs. (|35|); and 
we had seen how H could be written in terms of observables, Eq. (p8|). All we had to do was to refer back to those 
results. 

It should also be pointed out that this matrix formulation is very useful whenever we have non- uniform materials. 
For example, one might wish to study an experiment in which a kaon beam traverses vacuum, matter, and then vacuum 
again before it decays. Or a beam that traverses copper, carbon and then tungsten. In the matrix formulation, all 
we have to do is multiply three evolution matrices 

exp [—iHti] exp [—iH{t2 — ti)] exp [—iH{t3 — t2)], (55) 

each given by Eq. (|3^). 



IX. ON THE (MATHEMATICAL) RELATION WITH NEUTRINO OSCILLATIONS 

A. Boosted frames 

As we have mentioned before, the evolution equation (||) in which our study is based has been written in the rest 
frame of the P° — system. We denote this explicitly by 

. d 



dtr 



-l^iUest)) ^Hl^PiUest))- (56) 



The advantage of doing this is that in the rest frame the energy i s gi ven simply by E — m. As a result, the time 
parameter which appears in the solutions presented in Eqs. ( |3^ ) or (B7), through the time dependent functions g±(t) 
defined in Eq. (^), is really trost- 

Now imagine that we wished to have Eq. ( |56| ) given in a boosted frame (named the lab frame from now on) . In that 
case we would start by noticing that both the energy and the time are altered in the boosted frame. They become 

tlah = 7 ^rcst = — — ^rcst ■ (57) 

m 



Ignoring the matrix structure for the time being, Eq. (^6|) would change schematically into 

. d 



-lV^(ilab)) = m7 IV'(tlab)) = £^lab|V'(tlab)), (58) 



dtlah 

as it had to. Now, if the boost is much larger that the mass, p ^ m, we may use 



Euh ^^Jp^ + m?^p+ — + P+^ + ---- (59) 

However, we do not need to do this. We have already found the solution to Eq. ( |56| ) in the rest frame. In order 
to change it into the lab frame all we have to do is to substitute trost in the time evolution functions of Eq. ( p^ ) by 
^rest = i\ah/l- Wc noticc that 1/7 = m/E. Therefore, when written in terms of iiab the time evolution functions of 
Eq. (H) become 

5±(tr,,0 = e-^^^-Z^e-r-'-Z^^^) [ (go) 

i -isin(f ^ilab) 

And, using Am = ma — mb and m — {ma + mb)/2, we realize that the argument of the trigonometric functions is 
given by 

m A/i ml — ml i AT m 

-^—tuh = l^ilab- |— ^tlab. (61) 
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B. A neutrino-like oscillation 



For the comparison with neutrinos, it is most convenient to use the parametrization of the CP-violating quantities 
discussed in the first subsection of append ix To obtai n re lations that mimic those in the neutrino system, we 
compute the probability that P" becomes P° using Eqs. (13), d^), (il|), setting Im 4>ii = 1vd.9r = 0, and letting 
r = Ar ^ (another way of thinking about this limit is to suppose that we are performing an experiment in a time 
scale much smaller than the mesons' decay time). We find 



{P'\P\tu^)) 



sm Or sm 



(62) 



If, instead, the experiment is performed in matter, we obtain 



(P°|P"(tlab)) 



sm' 



2 a' 



sm 



hah , 



sin^ 9r 



|1 - 4r cos9fi + 4r2| 



4E 



sin Or 



(cos6'ij-2r)2+sin^ OrI 



where we have used Eq. 



in getting to the second line, and 



Ax Ax 



Ax/7 



E 



AE 



Ax 



2AA 2Am 2Am/7 ml 
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(63) 



(64) 



Eq. ( |63| ) exhibits a resonance structure because the time independent coefficient reaches its maximum if 2r = cos Or. 
For the final step in the connection to neutrinos, we look at this case further by assuming that the imaginary part of 
r (Ax), which is proportional to CTtot in appendix is negligible. Then Ax is real and we may parametrize 

ReAx 



V 



As a result. 



7 



EV 



is real and the resonance condition, which becomes 

2EV 



cos Or, 



(65) 



(66) 



(67) 



can be satisfied for 6r real. Eq. ( |63|) and the resonance condition in Eq. (^7|) are in exactly the same form as the 
usual discussions of neutrino oscillations in matter p8[ | . 

Although there is this mathematical connection between neutrino oscillations and pO — po oscillations, the situations 
are physically very different. Indeed, it is important to notice that there are no CPT relations between the two neutrino 
species involved in neutrino oscillation, and the vacuum mixing angle 6r in Eq. (^2|) may be small. Eq. (^) shows 
that, even if Or is small, the effective mixing angle in matter will be large when one hits the resonance condition in 
Eq. (|67|). In contrast, as we show in appendix in the P° — P'^ system, the deviation of cos Or from zero measures 
violations of CPT. Assuming CPT conservation in the P° — P° system, sin6'_R = 1 and the vacuum transition in 
Eq. ( p2| ) already reaches unity (at select times)^ Said otherwise, the small mixing angle mixing discussed in neutrino 
oscillations in vacuum, is (in the connection presented here) the mathematical analogue of large violations of CPT in 
the P" - po system. 



Recall that we have assumed V- 
multiplied by exp (— Ftrcst). 



=0 and Im(/)n = (T conservation). When F 7^ 0, the right hand side of Eq. 



appears 
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X. CONCLUSIONS 



We have shown that the presence of T violation in the neutral meson systems implies that the corresponding effective 
Hamiltonian H does not commute with itself. Therefore, H cannot be diagonalized by an unitary transformation 
and we must introduce the reciprocal basis. This basis must be used in order to obtain the correct form for some 
physical observables, such as the parameter A in the decays Bd — > J /ipKs- But, working with the reciprocal basis is 
a blessing rather than a nuisance. We show that using the reciprocal basis has the following advantages: 

• the relation between the effective Hamiltonian matrix when written in the mass and flavor basis is simply 
obtained and easily inverted, thus providing a parametrization of H in terms of measurable quantities; 

• one obtains a one line derivation of the evolution of the states; 

• propagation in matter is reduced to the case of propagation in vacuum, with the vacuum and matter parameters 
related in a trivial fashion, without any recourse to the Good equations; 

• the propagation in non-uniform media is reduced to a multiplication of evolution matrices. 

It is true that some of these results can be obtained without using the reciprocal basis. But, as we have tried to 
illustrate in the article, this concept is not only needed but, when used, greatly simplifies the various derivations. In 
addition, we can use this formalism to highlight the similarity between the matter effects in the — P° systems and 
the matter effects in neutrino oscillations. 
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APPENDIX A: MATTER EFFECTS IN THE REST AND LABORATORY FRAMES 



In this appendix we show how .ffnuc is related to physical cross-sections and what is the form of this relation in 
the rest and laboratory frames. This is relevant for Eq. (p9|) and for those wishing to expand on the analogy with the 



neutrino oscillations discussed at the end of section IX. We follow here the notation of references [E2l and E 



Let us consider the evolution of a coherent wave-packet with wavenumber k in the laboratory frame, 

-^(j) = ik(t). (Al) 
dz 

In the presence of a block of material at rest in the laboratory frame, the wavenumber suffers a shift given approxi- 
mately by 

fc'-fc-^/(0) = iv(^Re/(0) + z^), (A2) 

where N is the density of scattering centers in the medium and /(O) is the elastic forward scattering amplitude. On 
the last equality, we have used the fact that the imaginary part of /(O) is related to the total cross section atot by the 
optical theorem, 

Im/(0) = -^atot- (A3) 



We also recall that 



Ztt a cos tl 
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In this equation (and only here), 6 refers to the scattering angle in the laboratory frame. 
We conclude that the presence of matter changes the evolution in vacuum by an amount 

^-^c|>^-v{k'~k)^^~^f{0)<|>, (A5) 
diiab k/ V 

where v is the beam velocity in the lab frame and z — utiab- To change into the rest frame of the beam we notice 
that iiab = 7^icst and k — rwyv, where 7 — l/y/k'^ + m?, leading to 

^-^4>-~''-^m<t>■ (A6) 

direst m 

When studying the — P'^ systems, we denote by / (/) the elastic forward scattering amplitude of P° {P^)- 
Therefore, the new contribution in the — P'^ rest frame is given by 



(A7) 



where [[L7| 



leading to Eq. (E 



X = / and x = /, (A8) 



APPENDIX B: OTHER PARAMETRIZATIONS FOR T AND CPT VIOLATION 

The way we parametrize T and CPT violation in the mixing of neutral mesons is different from the parametrizations 
used by some other authors. For ease of reference, we collect here formulae summarizing the relationships among 
different parametrizations. 

1. The parameters (j)R and Or 

Some authors (for instance [l9|] ) introduce two complex angles 6ji and (f>ji by writing 

= 7V,cos^, g„ = TV^e'-^-sin^, 

= ArbSin^, = iVbe*-*«cos^. (Bl) 

Then, 

P 

6 ~ tanh (Im c/jr), 

e^-cosOR, (B2) 

and \/l — = sin6'ij. CPT is violated if and only if cosOr 7^ 0. T is violated if and only if Im^i^ 7^ 0. Some authors 
use a particular phase convention and claim that Re (j)R ^ also corresponds to T violation. Clearly this statement 
is false since the phase of q/p has no physical meaning; we know that there is one and only one T- and CP-violating 
quantity in H. 



With this notation, Eqs. (^, (|20|), (|2§), (g3|), and (|3^) become 

/ cos ^ - sin ^ \ / A'^. \ 

(B3) 



e''^«sin%- e'-^^cos^- / V ^ 
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(B4) 



H 



cos Or ^ 



e*"^« sin Or 



An 



- COS 9r^ 



(B5) 



( n,(t)+coseRq(t) e'"f"^ sine Rq_{t) \ / lpO\\ 

and 

|P°(t)) = [5+(t)+cos0flg_(i)]|P°)+e'*« sin0flg_(t) |/^), 

|P^(t)> - e-'^^- smf?;^..g_(t)|P") + [g+{t) - cosORg^it)] |P^), (B7) 

respectively. 

Finally, the relation between the matter and vacuum parameters described in Eqs. (js^) and (E^ become 





, X + X 






A/i' = 


V(AAt)2-2cos6ii? A^Ax4 


- (Ax)2 


= A^ - 4r cos Or + 4r2 


cos 9r = 


A/i cos 9r — Ax 




cos — 2r 


v/(A/i)2 - 2 cos 6*^ A^ Ax - 


t- (Ax)2 


VI - 4r cos^j^ + 4r2 ' 


sin0^ = 


Afj. suiOr 




sin (Jfl 


^{AfjiY -2 cos Or A^Ax- 


1- (Ax)2 


Vl - 4r cos^fl + 4r2 ' 



(B8) 

a/ — z COS £i/Li £ix + y^Xr V i — 'i''' COS UR i- ^r- 

and = (j)R. 



2. The parameters es and 5s 

Other authors (for instance |20y2l| ]) use two complex parameters, es and 5$, and write 

1 - es + (5s 



Obviously then. 



Pa 

qb_ 

Pb 



1 + es 



^5^ 
'^s 



1 + es + (^s 



(B9) 



q 
p 



^2 ■ 



8Re[ej(l + e|-<5|)] 



2<5s 



^1 



(l-e5)'-^l|) 



^1 



^2 ■ 



(BIO) 



CPT invariance corresponds to Js = 0. T invariance corresponds to Re [^s + 
this parametrization, however, always do so in conjunction with the assumption that Ss and es are small. Then, 



6g)] — 0. The authors who use 



2Rees, 
2ds. 



(Bll) 
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Moreover, y/1 - 9^ « 1 - 2(5|. 

It should be kept in mind that the i?-parametrization is exact and general, while the 5'-parametrization is interesting 
only when using a phase convention CV\P'^) = ±|P°), which implies that CP conservation corresponds to vanishing 
Ss and eg. 
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